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ABSTRACT 

We present and analyse near-infrared spectroscopy for a sample of 28 gravitationally- 
lensed star- forming galaxies in the redshift range 1.5 < z < 5, observed mostly with 
the Keck II telescope. With typical magnifications of ~1.5-4 magnitudes, our survey 
provides a valuable census of star formation rates, gas-phase metallicities and dynam- 
ical masses for a representative sample of low luminosity galaxies seen at a formative 
period in cosmic history. We find less evolution in the mass-metallicity relation com- 
pared to earlier work that focused on more luminous systems with z ~ 2 — 3, especially 
in the low mass (~ 10 9 M ) where our sample is ~ 0.25 dex more metal-rich. We in- 
terpret this offset as a result of the lower star formation rates (typically a factor of 
~10 lower) for a given stellar mass in our sub-luminous systems. Taking this effect 
into account, we conclude our objects are consistent with a fundamental metallicity 
relation recently proposed from unlensed observations. 

Key words: galaxies: evolution - galaxies: high redshift - galaxies: abundances - 
galaxies: kinematics and dynamics - gravitational lensing: strong 



1 INTRODUCTION 

The period corresponding to the redshift range 2 < z < 4 
is a formative one in the history of star- forming galaxies. 
During this era, mass assembly proceeds at its fastest rate 
and the bulk of the metals that establish present-day trends 
are likely manufactured. Through comprehensive multi- 
wavelength surveys, much has been lear ned about the demo- 
graphics of galaxies during this period (IShapley et all [2003 ; 
van Dokkum et al. l l2003l : IChapman et al. Il2005h as summa- 
rized in recent global measures (|Hopk ins fe Beacoml [20061 : 
lEllisll2008h . 

Attention is now focusing on the detailed properties of 
selected star- forming sources in this redshift range. Integral 
field unit (IFU) spectrographs aided by adaptive optics cor- 
rection on ground-based telescopes have delivered resolved 
velocity fields for various z > 2 galaxy samples revealing 



systemic rotation in a significant su b set (|Genzel et aD 
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spectroscopy, sampling rest-frame optical nebular emission 
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lines, defines the ongoing star formation rate and the gas - 
phase metallicity (jErb et al.1 [2006101 : iMannucci et alJl20Q9h . 
Metallicity has emerged as a key parameter since it mea- 
sures the fraction of baryonic material already converted 
into stars. Quantitative measures can thus be used to test 
feedback processes proposed to regulate star formation dur- 
ing an important period in cosmic history. 

A fundamental relation underpinning such studies 
is the mass-metallic ity relation first noted locally by 
iLequeux et al.l (Il979h and recen tly quantified in the SDSS 
survey by iTremonti et al ] (|2004h . The oft-quoted explana- 
tion for the relation invokes star- format ion driven outflows, 
e.g. from energetic super novae, which have a larger effect in 
low mass galaxies with weaker gravitational potentials. How- 
ever, other effects may enter, particularly at high redshifts 
where star formation timescales and feedback processes and 
their mass dependence likely differ. 

Motivated by the above, much observational effort 
has been invested to measure evolution in the mass- 
metallicity relation with redshift. The relationship has 
been defined using galaxy samples extending beyond z c± 
1 (lLamareille et al.l 1 20091: IPerez-Montero et al.1 l2009h and 



2 (|Erb et al.ll2006bl : lHallidav et al.ll2008l : lHavashi et al.1 
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l2009h . Most recently. IMannucci etHl (|2QQ9h have studied 
the properties of 10 Lyman-break galaxies (LBGs) at z ~ 3. 
These pioneering surveys have demonstrated clear evolution 
with metallicities that decrease at earlier times for a fixed 
stellar mass. 

Inevitably as one probes to higher redshift, it becomes 
progressively harder to maintain a useful dynamic range in 
the stellar mass and galaxy luminosity. In the case of most 
distant studies (e.g. IMannucci et al.ll2009h , only with long 
integrations can the mass-metallicity relation be extended 
down to masses of 1O 9 M , Samples defined via searches 
through gravitational lensing clusters are a much more ef- 
ficient probe of this important low mass regime. LBGs 
lensed by massive foreground clusters can be magnified 
by 2-3 magnitudes thereby probing intrinsically less mas- 
sive systems. Initial results using this technique have been 
presen t ed for small s amples by iLemoine-Busserqlle et al.l 
(|2003h . lHainline et al.1 (|20Q9h and lBian et al.1 (|2010h . " 

As part of a long-term program to determine the re- 
solved dynamical properties of sub-luminous z > 2 galaxies, 
we identified a large sample (~30) of gra vitationally- lensed 
systems with z >l-5 in the HST archi ve ([Sand et al.l 12005: 
ISmith et al.ll2005l : iRichard et alJlioiObh and embarked upon 
a systematic spectroscopic survey with the Keck II Near In- 
frared Spectrograph (NIRSPEC) to determine their emis- 
sion line characteristics. The initial motivation was to use 
the high efficiency of NIRSPEC to screen each target prior 
to more detailed follow-up with IFU spectrographs sampling 
the star formation r ate and velocity field across each source 
(| Jones et al.ll2010bL Livermore et al. in preparation). How- 
ever, a further product of this extensive spectroscopic sur- 
vey is detailed information on the star- format ion rate, emis- 
sion line ratios, and line widths for a large sample of lensed 
z > 1.5 galaxies. Our eventual sample comprises 28 objects 
including 5 from the literature mentioned above. The goal 
of this paper is thus to utilize this sample to extend studies 
of the mass-metallicity relation and related issues to more 
representative lower luminosity galaxies at early times. 

The paper is structured as follows. Section 2 introduces 
our sample and discusses the various NIRSPEC observations 
and their reductions as well as associated Spitzer data nec- 
essary to derive stellar masses. Section 3 discusses the mass- 
metallicity relation and the relationship between dynamical 
mass and stellar mass noting th at a subset of our sample 
has more detailed resolved data (| Jones et al.ll2010 b,a). We 
discuss the implications of our results in the context of mea- 
surements made of more luminous systems in Section 4. 

Throughout the paper, we assume a A-CDM cosmology 
with Qa=0.7, Q m =0.3 and h = 0.7. For this cosmology and 
at the typical redshift z ~ 2.5 of our sources, 1" on sky 
corresponds to ~8.2 kpc. All magnitudes are given in the 
AB system. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Lensed Sample 

We selected a sample of lensed galaxies for NIRSPEC follow- 
up using criteria similar to those used for fo r near- infrared 
spectroscopy of LBGs (e.g. lErb etaj]l2006bh but extended 
to lower intrinsic luminosities after correction for the lensing 
magnification. The relevant criteria are: 



• Availability of optical data from the Hubble Space Tele- 
scope indicating a prominent rest-frame UV continuum with 
V < 24 

• Spectroscopic redshift z > 1.5 derived from the lit- 
erature or as part of ou r Keck spectroscopic campaign 
([Richard et al.1 I2007L I2009L l2010bL Richard et al. 2010c in 
preparation) 

• An areal magnification factor \i > 1.5 mag provided by 
the foreground lensing cluster for which a well-constrained 
mass model enables a good und erstanding of the associated 
errors ( e.g. IRichard et aTll2010bl . Appendix A) 

• Emission lines predicted to lie in an uncontaminated 
region of the near-infrared night sky spectrum 

The application of these criteria generated a list of ~ 
50 arcs for further follow-up. We summarize in Table CD the 
23 sources drawn from this master list for which we were 
able to measure significant emission line fluxes. We have 
augmented this sample with ISAAC archival data and data 
from the literature for 5 other targets (see Sect. I2.2.2[) . In 
total, we consider data for 28 lensed sources spanning the 
redshift range 1.5< z <4.86. 



2.2 Near- infrared spectroscopic data 

2.2.1 NIRSPEC observations and data reduction 

The bulk of the spectrosc opic survey was condu cted with the 
NIRSPEC spectrograph ([McLean et alJll998h on the Keck 
II telescope in its low resolution mode during 7 observing 
runs (Table [2]). A 42 x 0.76" long slit was oriented along the 
major axis of each object, usually the direction of the highest 
magnification, in order to maximize the line fluxes (Fig. [T]). 
At this resolution, a different wavelength setup was selected 
for each target (filters Nl to N7, corresponding to the z 1 to 
K bands) for each group of lines ([Oil], H/3+[Oiii]aa4959,5007, 
Ha+[Nn] + [Sn]). 

A major advantage of targeting lensed sources is the 
ability to survey a large sample of low luminosity sources 
in an economic amount of observing time. We typically un- 
dertook 2 to 4 dithered exposures of 300 - 600 sees each, 
depending on the magnitude of the source and the sky lev- 
els in a given band. We used a three point dithering pattern 
with offsets larger than the size of the object along the slit. 
Standard stars were used as flux calibrators. Simultaneously 
with the NIRSPEC integrations, we took a series of 4-8 short 
exposures using the slit viewing camera, SCAM, in order to 
monitor the seeing and slit alignment with the object. 

The data was reduced using ID L scripts following the 
procedure described in more detail in Sta rk et al.l (J2007|) and 
IRichard et al.l (|2008h . Although the 2D spectrum is distorted 
on the detector, sky subtraction, wavelength and flux cali- 
bration were accomplished in the distorted frame, thereby 
mitig ating any deleterious effect of resampling (see iKelsonl 
2003 for more details). 

In com p arison with the general procedure applied by 
IStark et al.1 <|2007h for point sources, we took special care 
to prevent sky over-subtraction for our bright and extended 
sources. 
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Figure 1. Thumbnail HST images (V band) for each source targeted for near-infrared spectroscopy, showing the orientation of the 
corresponding NIRSPEC long slit. Targets are marked where there is otherwise ambiguity. 



2.2.2 Archival, literature and IFU data 

A modest amount of additional data on lensed galaxies is 
available in the literature, in the archive, or through new 
IFU data. To date, the Keck survey described above repre- 
sents the major advance. 

We included data from the ISAAC instrument 

on the ESQ VLT on 2 lensed galaxies at z—1.9 

(|Lemoine-Busserolle et all l2003h . Additional NIRSPEC 
data has been pub lished for two suitable objects by 
lHainline etHI (|2009h . Finally, we retrieved ISAAC archival 
data for the giant arc in C12244 ([Hammer et al.l [l989h at 
z — 2.24 which was studied by Lemoine-Busserolle et al. 
(2004). These archival J, H and K-band spectra have been 
reduce d with standard IRAF scripts, following the procedure 
used in lRichard et al. (2003). These 5 additional sources are 
listed separately in Table [1] 



To gether with the IFU data presented in I Jones et al.l 
(|2010bh , three of the targets described above (the cosmic 
eye, A1835 and MACS0451) were observed with the SIN- 
FONI integral field spectrograph on the VLT between May 
2009 and September 2010 as part of program 083.B-0108. 
In all three cases we used a 8x8" configuration at a spatial 
resolution of 0.25" /pixel and used J, H and if -band grat- 
ings which result in a spectral resolution of A/AA=4000 We 
used ABBA chop sequences while keeping the object inside 
the IFU at all times. Typical integration times were 7.2 ks 
(split into 600 second exposures) in <0.6" seeing and photo- 
metric conditions. Individual exposures were reduced using 
the SINFONI esorex data reduction pipeline and custom 
idl routines which, together, extracts, natfields, wavelength 
calibrates the data and forms the data cube. The final data 
cube was generated by aligning the individual data-cubes 
and then combining the using an average with a 3cr clip to 
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reject cosmic rays. For flux calibration, standard stars were 
observed each night during either immediately before or af- 
ter the science exposures. These were reduced in an identical 
manner to the science observations. Detailed analysis of the 
spatially resolved properties will be discussed in a forthcom- 
ing paper (Livermore et al. 2011 in preparation), but here 
we concentrate on the galaxy integrated emission line prop- 
erties measured from these observations. 
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2.2.3 Line fluxes and line widths 

Typical reduced spectra are shown in Figure [2] Line fluxes 
were measured using the IRAF task splot which uses both 
the science spectrum and the la error spectrum obtained 
from the extraction to derive the total flux and a bootstrap 
error. The line widths were measured on the highest signal- 
to-noise (> 10) spectra, and corrected for the effects of in- 
strumental resolution. Results are given in Table 2) We fit a 
single Gaussian to the Ha and [Oiii]A5007 lines, and a dou- 
ble Gaussian to [On] with the lines fixed at rest wavelengths 
of 3726.1,3728.8 A. We constrain the [On] lines to have the 
same width and an intensity ratio 7(3726. l)//(3728. 8) ~ 1 
as seen in high red shift galaxies observed with higher spec- 
tral resolution (e.g. ISwinbank et al.l l2009). 

The spectral resolution R — X/FWHM measured from 
bright, unblended OH sky lines in the NIRSPEC spectra is 
found to vary linearly with spectral order m as expected, 
and the ratio R/m varies smoothly with wavelength from 
roughly 2100 at A = 1 fun (m = 4) to 4000 at 2.2 /mi (m 
2). This results in an instrumental FWHM ranging from 
160-260 km/s (a = 70-110 km/s) for the lines observed 
in this work. Measured line widths exceed the instrumental 
resolution in all but one case (RXJ1347-11), for which the 
lcr upper bound is given in Table 2) 

The uncertainty in the line widths is propagated from 
the la error spectrum. Gaussian fits used to determine the 
line width have residuals xl — 1 suggesting that the er- 
ror spectrum is a good estimate of the intrinsic noise. After 
subtracting the instrumental resolution in quadrature, line 
widths are typically determined to ~10% accuracy. The un- 
certainty is significantly larger in cases where the line of in- 
terest is blended with a sky line, or has a width close to the 
instrumental resolution. We estimate any systematic uncer- 
tainty in the measured instrumental resolution to be < 3% 
based on independent measurements of multiple bright sky 
lines, hence we expect measurement errors to dominate. 

We can test whether the velocity dispersions mea- 
sured with NIRSPEC are reliable, where there is over- 
lap with integral field data, by undertaking a comparison 
with the more extensive 2D velocity data (ILaw et al.l [2009 ; 
iForster Schreiber et all 120091 : Ijones et al.l l2010bh . Here we 
enlarge the comp arison sample b y considering all relevant 
NIRSPEC data (jErb et al.ll20Q6bh . As Fig.[3]shows, the data 
is in general agreement with the mean IFU velocity disper- 
sion being 1.2 ± 0.3 times that measured with NIRSPEC. 
This suggests that the longslit data provide a reasonable 
estimate of the global galaxy dynamics and cab be used, 
for example, with the spatial extent to measure dynamical 
masses. 
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Figure 3. Comparison of the velocity dispersion measured from 
our long slit NIRSPEC data with 2D integral field unit data from 
various sources. Blue points refer to mean d ispersions (a) from 
OSIRIS IFU observations reported in ljones et al.l (l2010blV C reen 
crosses show equivalent OSIRIS data from lLaw et al . (2009) and 
red c rosses are integrated velocity dispersions from the SINFONI 
IFU (Fors ter Schreiber et ail feoOQV NIRSPEC measurements for 
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the latter samples are reported in lErb et al . (2006b). 



2.2.4 Line ratios from stacked spectra 

One of our aims is to measure the mass-metallicity relation 
in our sample, extending it to low- luminosity sources. Given 
the short exposure times for our survey which has enabled 
such a large sample to be constructed, inevitably many of 
the fainter diagnostic lines, e.g. [Nil] and [Sll] are not always 
visible in individual spectra. In order to estimate the typ- 
ical prevalence of such faint emission lines, we construct a 
stacked spectrum about the Ha line. 

This has been done for two cases. Firstly, for all 9 lensed 
sources where Ha is detected at S/N > 10, and, secondly, 
for all 6 cases where the stellar mass (§2.4) is < 3x 10 9 M . 
First we resample the Ha spectra into the rest frame with 
a common dispersion of 1 A, and then scale each spectra to 
a common Ha flux. We reject the minimum and maximum 
values at each wavelength in order to remove outliers (due to 
sky line residuals, for example), and average the remaining 
data. The results are consistent with Gaussian noise. The 
composite spectra are shown in Figure |4] 

We measure line ratios in the stacked spectra by fitting 
a Gaussian profile to Ha, then using the position and line 
width to determine the [Nil] and [Sn] fluxes and their boot- 
strap error. The resulting ratios are provided at the bottom 
of Table E] 



2.2.5 Aperture corrections and further checks 

One obvious limitation of long-slit spectroscopy is a signifi- 
cant fraction of the line flux may be missing, especially for 
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Figure 2. Typical extracted NIRSPEC spectra (flux units in erg s 
key nebular lines were covered using 3 different spectrograph setting 
dotted line indicates the la error on the spectral flux. 
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Six sources in our sample have been presented by 
Ijones et al ] (|2010bh and were therefore covered both with 
NIRSPEC and OSIRIS. This allows us to perform a further 
check on the aperture correction factors. We found that the 
total line fluxes agree for each object within 10%, which is 
likely due to the aperture correction and/or the relative flux 
calibrations. We adopt this 10% error estimate on the total 
line flux when estimating the SFRs later in the paper. 

As mentioned earlier, we can use the same six targets 
to compare the velocity dispersions estimated from the NIR- 
SPEC and OSIRIS data on the same emission lines (§2.2.3). 
Although the long slit estimates have larger error bars, they 
are consistent (within la) with the more reliable IFU esti- 
mates. Nonetheless, we adopt the IFU value whenever avail- 
ble. 
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Figure 4. Stacked spectra chosen to reveal fainter emission lines. 
Top: stack of all lensed sources with S/N > 10 in the Ha line. 
Bottom: stack of lensed sources with stellar mass < 3 x 10 9 Mq. 
Dotted lines show the positions of Ha, [Nil] AA6548, 6584, and 
[Sn]AA6717, 6731. 



the most extended objects or in poor seeing conditions. Our 
slit width was 0.76 arcsec and so the latter may be a con- 
cern for some objects (see Table [2]). In order to compare 
measurements based on long-slit spectroscopy with photo- 
metric estimates derived from imaging (see Sect. 12. 4p . we 
carefully estimated aperture correction factors. 

The procedure adopted was the following: we used the 
SExtractor segmentation map to select those pixels of the 
HST image associated with the object. We smoothed this 
new image by our estimate of the seeing conditions (Table [2} 
and measured the fraction of the total flux falling inside the 
region covered by the long slit. We found aperture correction 
factors between 1.2 and 2.9 depending on the object. The er- 
ror on these corrections was estimated using a iO^'error on 
the location of the slit, as checked from the SCAM images. 
The correction factors are later used to derive the total star 
formation rate (SFR) of the object, assuming the measured 
equivalent width of the lines are representative of its average 
across the entire galaxy. Of course, the metallicity and line 
ratio measurements are unaffected as they rely entirely on 
the spectroscopic data. We also checked that measurements 
of the 8 o'clock arc are con sistent with the results presented 
bv lFinkelstein et al.1 (|2009h . 



2.3 Magnification and source reconstruction 

In order to correct all physical properties (star- format ion 
rates, masses, physical scales) for the lensing magnification, 
we correct each source for its corresponding magnification 
factor \i which stretches the physical scales of each object 
while keeping its surface brightness fixed. As discussed in 
§2.1, our targets were specifically chosen to lie in the fields 
of clusters whose mass models are well-constrained from as- 
sociated spectroscopy of many lensed sources and multiple- 
images (see example in Appendix A and I Richard et al.l 
(| 201 Obi ) for further discussion). 

The values of \i are obtained through modeling of the 
cluster mass distribution using the Lenstool software El (see 
references in Table [J), or from the literature (for the addi- 
tional targets in Table [1]). The relative error on /i, when de- 
rived from a Lenst ool parametric m odel, follows a Bayesian 
MCMC sampler (| Julio et al.ll2007l ) which analyses a family 
of models fitting the constraints on the multiple images. In 
Table [1] we report the final values of and their associated 
errors; these are used to correct /i-dependent physical pa- 
rameters (star-formation rates and masses) throughout the 
rest of the paper. Note that the magnification factor does 
not affect any of the line ratios measurements. 

We also construct the demagnified (unlensed) source 
morphology of each arc modeled with Lenstool by ray- 
tracing the high-resolution HST image back to the source 
plane, using the best fit lens model. By comparing the sizes 
of the observed and reconstructed image of a given target, 
we can verify the agreement of the m agnification factor used 
and this size ratio (c.f. discussion in I Jones et al.l [201 Obh . 



2.4 Photometric measurements 

A large variety of space- and ground-based images are avail- 
able for each object from archival sources. All HST im- 
ages (optical and near- infr ared) have been reduc ed using 
the multidrizzle package (Koe kemoer et al . 2002), as well 
as specific IRAF script s for NICMOS data, as described 
m iRichard et all (2008). Ground-based near-infrared im- 
ages were re duced following the f ull reduction procedure 
described in I Richard et al.1 (|2006l ) , and calibrated using 
2MASS stars identified in the field. 



1 |http: //www. oamp.fr/c osmology/lenstool/ 
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Archival IRAC data in the first 2 channels (3.6 and 4.5 
/zm) are available for all sources except MACS 1423, but we 
only consider sources which are not contaminated by nearby 
bright galaxies when deriving their photometry. We combine 
the post-BCD (Basic Calibrated Data) frames resampled to 
a pixel scale of 0.6". 

The HST image providing the highest signal-to-noise 
was used to measure the i ntegr ated brightness with SEx- 
tractor (Bert i~fe Arnoutsl Il996h , whereas the double im- 
age mode was used to measure the relative HST colours 
inside a l" aperture. A small aperture correction is applied 
to all HST photometry to deal with the PSF differences be- 
tween the different ACS, WFPC2 and NICMOS bands. In 
the case of ground-based colours, the primary HST image 
was smoothed by a Gaussian kernel corresponding to the 
measured seeing, and convolved with the IRAC PSF (de- 
rived from bright unsaturated stars in the image) in order 
to incorporate IRAC colours. The full photometry is sum- 
marised in Table lA2l in appendix. 



3 PHYSICAL PROPERTIES OF THE SAMPLE 

Table 2] summarises the derived quantities which will form 
the basis of our analysis. We now discuss the physical mea- 
sures in turn. 

3.1 Star formation rate and AGN contribution 

The intrinsic star formation rate (SFR) is estimated from 
the total (aperture-corrected) flux in the Balmer lines (fua 
and fnd) based on the well- constrained calibrations by 
IKennicuttI (|l998l) , including the correction for magnifica- 
tion factor fi and its associated error. In the absence of 
Ha we assume a typical ratio Ha/H/3=2.86. For two ob- 
jects (RXJ1347 and A1689-Highz), neither Ha nor R/3 is 
available, therefore we use the [Oil] line to derive the SFR, 
although it is a less robust estimator because of metallicity 
and excitation effects fe.g. iGallagher et al.lll989h . 

The derived values of SFR span a wide range, between 
0.4 and 50 M /yr, and are typically a factor of 10 lower than 
other samples of LBGs with emission line measurements 
(Fig. [5]). The majority (16 objects) of the lensed sources 
have S FR< 4M / yr, va lues whi ch are absent from t he sam- 
ples of lErb et all |2006bl ) and lMaiolino et al.1 (|2008l ). 

An independent estimate of the SFR can be obtained for 
every object using the ultraviolet continuum luminosity (L u ) 
estimated at 1500 A rest-frame. We measure L v by fitting a 
power-law fx oc A _/3 to the broad-band photometry between 
1500 and 4000 A rest-frame. The UV slope (f3) is given in 
the first column of Table IA2I A value f3 = 2.0 (constant 
AB magnitude) is assumed when only a single photometric 
datapoint was available. T he SFR derived from the UV- 
continuum is then given bv IKennicuttI (1998). 

We can determine from nebular line ratios whether 
there is a strong contribution from active galaxy nuclei 
(AGN) in the line emiss i ons. T his is done through the di- 
agram of lBaldwin et~all £981) (BPT dia gram), where we 
overplot the location of our targets in Fig. [6j 

The majority of our sources lie in the region of the BPT 
diagram where sta r- forming galaxies are commonly found 
(|Kewlev et al.ll200lh . MACS0451 and AC114-A2, however, 
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Figure 6. Location of the NIRSPEC targets (blue diamonds) 
over the BPT diagram, used as a diagnostic for AGN vs star- 
formation (see text for details). The dotted line is the em- 
pirical separation betw een star-forming and AGN objects from 
IKauffmann et al.l (120 03), and the dashed line the theoretical sep- 
aration from lKewlev et al l (l200lh . Other sources from the litera- 
ture are shown as green points (lense d objects fromlHamline et ahl 
120091 and lemoine03) or red triangles (|Erb et al.l2006b|j201oh . The 
two circled points are objects showing signs of AGN activity, as 
discussed within these papers. 

lie close to the boundary with the region occupied by AGNs. 
We consider it unlikely, however, that these objects are 
AGN-dominated, because of the lack of X-ray emission in 
Chandra data and Civ lines in the optical spectra. Fur- 
thermore, the Spitzer data shows evidence for the 1.6/im 
rest-frame stellar bump and no sign of obscured AGN activ- 
ity, which woul d produce a raising s lope in the redder IRAC 
channels fe.g.. Illainline et al.ll201oh . The location of these 
objects in the BPT plane simply suggests a higher radiation 
field, similar to that seen in loc a l starb ursts and the sample 
of LBGs studied by lErb et al.1 (|2006al ) at z ~ 2 (see also 
Sect 13. 5|) . We thus conclude that the nebular emission we 
see in our sources arises from intense star formation. 



3.2 Extinction 

Dust extinction plays an important role when deriving the 
physical properties of galaxies, as it will affect the observed 
line fluxes and some of the line ratios. One of the estima- 
tors we can use to measure this extinction is the UV spectral 
slope /3, which is related to the extinction affecting the young 
stars. We can also use the ratio between the two SFR esti- 
mates (from the observed UV continuum and the Ha emis- 
sion lines) as an independent estimator, as it reflects the 
differential extinction between the two wavelengths. Note 
that even in the case of a perfect reddening correction, the 
nebular emission and t he UV do not pr obe star- format ion 
on the same timescales (|Kennicutt| [l998) . 
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Figure 5. Probability distribut ion of star-form ation rates in our sample of lensed objects (blue histogram) co mpared to the sample s of 
LBGs at z ~ 2 (red histogram, lErb et alj[2006bl ) and at z ~ 3 (dotted histogram, AMAZE and LSD samples Mai olino et al"ll2008l and 
iMannucci et al. II 2009). The left panel compares the distribution of SFR before correcting for extinction, while the right panel compares 
the values after applying the extinction correction. 



Figure [71 compares both estimates for a subsample of 
the l ensed objects , toge ther with the relation predicted by 
the lCalzetti et all (|2000l ) extinction law. It shows that there 
is a general agreement with the theoretical predictions, al- 
though with quite a large scatter. One of the reasons for the 
differences is probably a different extinction factor affecting 
the young and old stellar populations, or a measurement bias 
towards low extinction regions when measuring /3. This is il- 
lustrated by the location of the sub-millimetre lensed galaxy 
A2218-smm (|Kneib et al.ll2004n . which has a high extinction 
but a measured slope f3 = 1.6 (white diamond in Fig. [Tj) . 

One way to overcome this issue is to use the full SED 
(from rest-frame UV to near-infrared) in orde r to derive the 
best e xtinction estimate E(B-V) assuming the lCalzetti et al] 
(2000) law (see Sect l3.3.T]) . 8 objects in our sample also have 
both H a and H/3 detected, and the corresponding values of 
the Balmer decrement range between 2.5 and 4.5, consistent 
with < E(B — V) < 0.35 although with large uncertainties. 

We use the best-fit E(B-V) values, given in Table H to 
correct the SFRs, individual line ratios and metallicities in 
the next sections. In particular, reddening has a strong ef- 
fect on individual SFRs, and it is very sensitive to unknown 
factors (differential reddening in a given object, presence 
of dust in the lightpath through the galaxy c luster, devi- 
ation from the Calzetti law as pointed out by ISiana et~al] 
2009). We present both the uncorrected and the extinction- 
corrected (SFR corr ) in Table [H and when comparing with 
other samples (Fig. [5]). 
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Figure 7. Comparison between the SFR estimated from the UV 
continuum and the Ha emission line, as a function of the UV slope 
p. The solid curve is the theoretical prediction for t he relation 
between the two values using the Calzetti et al . (2000) extinction 
law. 



measuring the widths of the most prominent spectral lines 
allow us to infer dynamical masses, which should be closer 
to the total baryonic mass of these galaxies. 



3.3 Masses 

Two mass estimators can be derived from the available data. 
Multi- wavelength broad-band photometry gives us access to 
the stellar mass, through the modelling of the SED, while 



3.3.1 Stellar masses 

Our stellar masses are derived using the precepts discussed 
in detail bv lStark et al.l (|2009h . We derive the stellar masses 
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for our sample by fitting the Chariot & Bruzual (2007, S. 
Chariot, private communication) stellar population synthe- 
sis models to the observed SEDs. We consider exponentially- 
decaying star formation histories with the form SFR(t) ~ 
exp(— tr) with e-folding times of r = 10, 70, 100, 300, and 
500 Myr in addition to models with continuous star forma- 
tion (CSF). For a given galaxy, we consider models rang- 
ing in age from 10 Myr t o the age of th e Universe at the 
galaxy's redshift. We use a | Salpetei (1 19551) initial mass func- 
tion (IMF) and the ICalzetti et al.1 (|2000h dust extinction 
law. Finally, we allow the metallicity to vary between solar 
(Z0) and 0.2 Z©, the range found for the gas-phase metallic- 
ity using nebular line ratios (Sect. [33]). We account for th e 
intergalactic medium absorption following iMeiksinl (|2006h . 
The best fit values of the stellar mass M* and extinction 
E(B-V) are summarised in Tabled 

We note that the presence of strong emission lines may 
affect the broad-band photometry and therefore the mass es- 
timates. We estimate the contribution of the strongest emis- 
sion lines affecting the H and K band magnitudes. In the 
most extreme cases (largest equivalent widths) the H and K 
band flux would be affected by at most 5-10%, which does 
not have a significant effect compared to our estimated er- 



3.3.2 Dynamical masses 

We now use the velocity dispersions a (measured in Sect. 
I2.2.3P to estimate of the dynamical mass. We will express 
the virial masses Md yn of our objects as a function of a and 
their typical size. 

Assuming the ideal ized case of a sphere of uniform den- 
sity (|Pettini et al.ll200lh we have: 

Mdyn = Cr 1 / 2 cr 2 /G with C — 5 or more conveniently, 



M d y 



1.16 x 10 1U Mrr 



n/2 



(100 km s- 1 ) 2 kpc 



(1) 



where 7*1/2 is the half-light radius, which we measure on 
the source plane reconstructions of our targets (see Sect. 12. 3[) 
using the FLUX_RADIUS parameter from SExtractor. This pa- 
rameter estimates a circularized size corresponding to half of 
the total detected fluxes. Thanks to the high magnification 
our sources are well resolved in the HST images, but the 
main source of error in estimating 7*1/2 (which has a strong 
impact on Mdyn) is the source reconstruction itself. We pro- 
duced 100 reconstructions of each source sampling the dif- 
ferent parameters of the lens model, using the MCMC sam- 
pler described in Sect. 12.31 and use them to derive the mean 
and dispersion on the measured 7*1/2 . We also checked this 
measurement independently using a half-light radius defined 
with the Petrosian radius at 20% o f the central flux of th e 
object, similar to the work done by ISwinbank et all (|2010h , 
and found consistent results. 

The values of 7*1/2 and Md yn are reported in Table[4] We 
note that the geometric correction factor C can vary signifi- 
cantly and is likel y between 3 a nd 10 for these highly turbu- 
lent galaxies (|Erb et al.ll2006bh . Hence the true dynamical 
mass can differ by up to a factor of 2, while uncertainty in 
the radius and a is much smaller. Our choice of C = 5 gives 
slightly highe r dynamical masses than in some othe r studies 
(e.g. C 3.4. lErb et al.ll2006bl : iBouche et alJl2007h . 



3. 3. 3 Comparison 

The range of stellar masses spanned by our lensed objects 
(typically 5 10 8 - 5 IP 1 " M^ is l o wer by a factor of 3-5 
than those surveyed by Erb et al.1 ([2006a ) ; iMaiolino et al.l 
(|2008T ) and lMannucci et al.1 (| 200$ ) ( Fie:. 151 left). Within the 
subsample with reliable dynamical masses, we can directly 
compare the dynamical and stellar masses. Although we dis- 
cover a large dispersion (Fig. [51 right), the average ratio 
< log(Md yn /M*) >= 0.36 dex is a factor of ~ 2, suggest- 
ing of ^40% of dynamical mass not prese nt into st a rs. Thi s 
value is similar to the result found by lErb et all (|2006d ) , 
and confirms the dominance of baryonic mass in the cen- 
tral regions of t hese g alaxies. This was already pointed out 
by IStark et al.1 (2008) in the case of the Cosmic Eye. We 
note that a different choice for the multiplicative factor in 
Eq.[T] would give ~ 30% lower dynamical masses and would 
strengthen this conclusion. 

The gas mass can be evaluated from the SFR and the 
size of the objects 7*1/2 (ideally assuming the star- format ion 
is unifo rm over the projec ted surface seen in the UV) , assum- 
ing the iKennicuttl (| 19981 ) law between star-formation rate 
and gas densities applies at these redshifts. We use this re- 
lation to derive the total gas mass: 

M gas = 5.03 10 8 SFR " 71 r?/ 5 2 8 M (2) 

Defining the gas fraction f s as aS fgas — -^fgas / ( -^gas ~~\~ 
M*) , we find gas fractions ranging from 0.2 to 0.6 after 
reddening correction, compatible with what is derived from 
the comparison of stellar and dynamical masses. This is 
in average ^50% lower that the gas fraction measured by 
lErb et al.l (|2006ch at z ~ 2 and bv lMannucci et al.1 (|2009h in 
the Lyman-break galaxies Stellar populations and Dynamics 
(LSD) sample at z ~ 3, and reflects the lower SFRs in our 
objects for a given stellar mass. However, we see the same 
trend of a lower / gas towards the higher masses. 



3.4 Metallicities 

The relationship between the oxygen abundance, or gas 
metallicity Z defined as log(Z) = 12+log(0/H), has been ac- 
curately calibrated against line ratio s of the prominent nebu- 
lar emission lines (Nag ao et al.l l2006). Suitable line ratios are 
available for half of the objects in our sample, and depend- 
ing on the availability we can combine the estimates from 
three different metallicity d iagnostics, usin g the r ecent em- 
pirical relations deriv ed by IMaiolino et aTl |2008) (see also 
iMannucci et alJl2009h . 

The prime oxygen abundance indicator is the R23 ra- 
tio (|Pagel et al.lll979h : R23=([On] + [Oin]AA4959,5007)/H/3. 
As the [Oiii]A4959 line is usually weakly detected we adopt 
a canonical value of 0.28 for the [Oiii]A4959/[Oin]A5007 
ratio. The second estimate for oxygen abundance is the 
032 ratio: O32=[Om]A5007/[Oii]. The final estimate is the 
N2=[Nn]/Ha indicator. 

Defining the parameter x = IjxA Z/Zq) and using 
/o#(Z©)=8.69 (jAllende Prieto et alJlioOll), w e adopt the fol- 
lowing fitting formula (Mai olino et al.ll2008T ) on the redden- 
ing corrected line ratios: 
log(R23) = 



0.7462 - 0.7149 x - 0.9401 x 2 - 0.6154 x 3 



0.2524 x 4 (3) 
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Figure 8. (Left) Probability di stribution of s tellar masses in our sample of lensed objects (blue histogram) compared to the samples of 
LBGs at z ~ 2 (red h istogram, lErb et al.ll2006ch and at z ~ 3 (dotted histogram, AMAZE and LSD samples, iMaiolino et al.ll2008l and 
iMannucci et al.ll2009l ). (Right) Comparison between stellar and dynamical mass estimates. In general we find a higher dynamical mass 
compared to the stellar mass, suggesting the presence of large amounts of gas mass in this sample. 



log(032) = -0.2839 - 1.3881 x - 0.3172 x 2 



(4) 



log(N2)= 

- 0.7732 + 1.2357 x - 0.2811 x 2 - 0.7201 x 3 - 0.3330 x 4 (5) 

These equations provide three different estimates for Z: 
Zr23, ^032 and Zm- Of these, the Z N 2 and Zr23 estimates are 
the ones showing the lowest dispersion (typically 0.1 dex) 
while the Z032 estimate has a dispersion of 0.2-0.3 dex. How- 
ever, the Zr23 relat ion has a "two-branch degeneracy" (e.g. 
IPettini et alJ l200ll ) between a low-metallicity and a high- 
met allicity value. In order to discriminate between the two 
values, we use the Z N2 calibration when available (usually 
synonym of a source on the upper- metallicity branch), oth- 
erwise the Zq32 calibration. For the 5 sources taken from the 
literature, we recalculate the best fit Z using the published 
line fluxes and reddening factors. The resulting value of Z 
is consiste nt with the published value exce pt in the case of 
AC114-S2 lLemoine-Busserolle et all <|2003h , where the pub- 
lished value log(Z) — 7.25 ± 0.2 makes this object a clear 
outlier in the mass-metallicity relation (Sect. I4.ip . Taking 
the published line fluxes, the N2 limit derived places it at 
the junction between the two branches of the R23 relation. 
To resolve this ambiguity we use the values derived with the 
upper-branch. 

The final best- fit value for log(Z) is reported in Tabled 
for each source. The error bars include the typical dispersion 
in the relation. We find metallicity values ranging from 0.25 
Z to 1.7 Z . 



3.5 Ionization parameter 

In the previous section, the 032 parameter has been used 
mainly to discriminate between the lower and upper-branch 
of the R23 calibration of oxygen abundance. For a given 



metallicity, the 032 line ratio can also be compared with 
models of Hn regions to measure the ionization parameter 
U, i.e. the ratio of densi ty of ionizing photons ov er the den- 
sity of hydrogen atoms ((Kewlev &; D opita 20021 ). Bas ed on 
early samples of bright LBGs from lPettini et all (l2Q0lh , high 
i oniza tion levels have been measured by iBrinchmann et al.l 
(2008), with \og(U) ~ —2.0, compared to local samples. This 
same ionization parameter shifts the objects upward in the 
BPT diagram (Fig. [6}, perhaps p ushing the frontier between 
star- forming galaxies and AGNs (|Erb et al.l l2006b). 

The 032 parameter can be measured for 6 objects in our 
sample, as well as 2 upper limits, and we obtain values of 
032 in the range <032< 1. For the range of metallicities 
foun d previously (0.2 to 2.0 Z©), we use the curves provided 
bv lKewlev fe Dopital {2002) to derive ionization parameters 
— 2.9 < log (17) < — 2 whereas the typical values from local 
galax ies are in the range — 4 < log (£7) < —3 ([Lilly et al.l 
l2003h . 

This can be illustrated by constructing the 032 vs R23 
diagram, as the R23 metallicity estimator has only a weak 
dependance on t he ionization p arameter. We overplot the 
results found by lHainline et al.l (|2009h on 4 lensed objects 
(included in our sample) and add our 6 new measurements 
and 2 upper limits in this diagram (Fig.[9])- We can see that 
, in average, our sample is systematically shifted towards 
high er values of 032 co mpared to the lower redshift objects 
from iLillv et al.l ((2003). A likely explanation for this effect 
might be that the physical conditions in the relevant Hn 
regions are different from those in the local Universe, for 
example, with larger electron d ensity and/or larger escape 
fraction (|Brinchmann et al . 200i). 

An even more ext reme resul t was found recently in a 
high redshift object bv lErb et al.l (|2010h , where they derive 
an ionization parameter log{U) ~ 1.0. The very young age 
found in this object is one of the factors explaining such a 
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Figure 9. 032 vs R23 diagram, showing the effect of the ion - 
ization parameter on 032 (adapted from Hain line et all 2009). 
The red triang les are the low redshift (z < 1) sample of galax- 
ies from iLillv et al.l (120031 ). Our sample of z >1.5 lensed galaxies 
is shown in bl u e dia monds, while other literature data used by 
Hainline et al . (2009) is shown in green. Sources in our sample 
showing young stellar populations (best age < 100 Myrs from the 
SED fitting) are circled and typically lie in the top part of the 
diagram. 

high value of U. Indeed, we can see some trend with age in 
the 032 vs R23 diagram, despite the small number of objects 
in our sample. By selecting sources having very young stellar 
populations (best age < 100 Myrs from our SED fitting), 
they all lie in the top part of this diagram, with the highest 
032 values. 

4 THE MASS METALLICITY RELATION 

4.1 Comparison with Earlier Work 

We now compare our measurement of the mass-metallicity 
relation at high redshift with that of earlier workers recogniz- 
ing that our data, for the first time, probes to lower masses 
and lower SFRs by virtue of our selection of gravitationally- 
lensed systems (Fig. [5] and [8]). Figure [10] summarizes the 
current situation. We find that all galaxies in our sample 
lie below the well-den ned metallicity relation at z — 0.07 
(|Tremonti et alJ l2004|) therefore supporti ng strongly the 
case for evolution. IMannucci et all (|2009l ) have proposed 
two best fits for th i s relat ion, at z = 2.20 and z — 3.0, 
from the lErb et al.l (l2006ah. the AMAZE (jMaiolino et all 
l2Q08h and the LSD jMannucci et al.1 l2009h samples. If we 
partition our sample into two redshift bins: 1.5 < z < 2.5 
and 2.5 < z < 3.5, we can see that, on average, the mass- 
metallicity relation follows these trends. However we observe 
a large scatter (~ 0.25 dex) in the 1.5 < z < 2.5 range and 
our the evolution implied by our d ata is less extreme ove r 
2.20< z <3.3 than that suggested bv lMannucci et al.l (|2QQ9h . 



Compared with the earlier unlensed studies, our sample 
spans a wider range of stellar masses. Therefore, even if our 
median stellar mass i s similar to the z ~ 3 sample from 
IMannucci et al.l {[2009) , we have significantly increased the 
number of galaxies at M ~ 10 9 Mq for the most highly 
magnified objects. At these masses, we do not seem to see a 
steep decline in metallicity seen in the z = 2.2 and z = 3 best 
fit found in the earlier surveys (Figure HP)) . Instead, we see 
an average shift of ^0.25 dex towards higher metallicities for 
low-mass objects (~ 10 9 M©). Although this trend is limited 
by the small statistics of our sample, it is also visible when 
overplotting the results from the composite spectra denned 
inSect.EOO 

4.2 A fundamental metallicity relation? 

By comparing the histograms in Fig. [5] and Fig. [8] it is clear 
that our lensed objects have much lower SFRs than those in 
more luminous LBGs, even when accounting for the redden- 
ing correction, but only slightly smaller masses. Indeed, the 
effect of the increasing SFR at higher redshift is a usual ex- 
planation for the evolution of the mass-metallicity relation 
seen in Fig. 1101 and selecting objects of lower SFR would 
remove this effect and explain the slight deviation of our 
sample towards slightly higher metallicity. The reason for 
the much lower SFR in our sample is expected, as the un- 
lensed samples of LBGs are selected through their SFR in 
the UV, while we made no strong assumptions on this pa- 
rameter. 

Recently, IMannucci et al.1 (2010) have proposed to in- 
clude the SFR as a third component in the mass-metallicity 
relation, which would explain its evolution in redshift. One of 
the reason for the influence of star- format ion is that outflows 
would be more efficient in low-mass objects. Exploring the 
third-dimensional space denned by stellar mass, gas metal- 
licity and SFR, they fit a surface in low redshift galaxies 
denned by: 

Z est = 8.90+0.37 ra-0.14 s-0.19 ra 2 +0.12 m s-0.054 s 2 (6) 

with m log(M*) - 10 and s log (SFR). 

By comparing the measured metallicity Z with the pro- 
jection Zest estimated from this surface, they managed to 
reduce the dispersion in SDSS galaxies to ~ 0.05 dex. Fur- 
ther computing the distance of higher redshift samples from 
this surface, they found no significant evolution in this fun- 
damental relation upto z ~ 2.5, within the la error. 

We computed the values of Z es t for each galaxy in our 
sample and compare it with the measured Z. We find that 
in average, the best-fit surface predicts the metallicity with 
no significant offset, and with a scatter of ~ 0.2 dex, making 
this metallicity relation compatible with our measurements 
within a 2a level. When plotting the distance to the fitted 
surface as a function of stellar mass (Figure 1 1 If) we see no 
significant trend with stellar mass in both redshift ranges 
probed by our sample. 

4.3 Summary and perspectives 

We have presented the results of a near-infrared spec- 
troscopy survey targetting 23 bright lensed galaxies at z > 
1.5, complemented by 5 sources from the literature. We sum- 
marize here our findings: 
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Fi gure 10. The relation between gas metallicity and stellar mass as a function of redshift. This is based on the diagram presented 
bv lMannucci et all ((2009), with their estimates on z ~ 3 galaxies (LSD sample) shown as grey datapoints. The local mass-metallicity 
relation is showed as a dotted line for z = 0.07, while their best fit at z ~ 2.2 and z ~ 3.0 (including both the LSD and AMAZE samples) 
are shown as solid lines in the left and right panels respectively. We present as diamond symbols the subsample of 1.5 < z < 2.5 lensed 
galaxies in the left panel, and the 2.5 < z < 3.5 sample in the right panel. The values derived from the two co mposite spectra (see Sect. 
12. 2.41) are shown as red points in the left panel. The original value of 12 + log(0/H) for AC114-S2 published bv lLemoine- Busseroll e et al.l 
(|2003h is given as a triangle. 



• After correction for the magnification factor, our sam- 
ple shows in average 10 x smaller star- format ion rates and 
5x smaller stellar masses than the samples of LBGs at the 
same redshifts. Such low values of SFR would not be acces- 
sible without the strong lensing effect, making our sample 
complementary to LBG studies. 

• The comparison of dynamical and stellar mass esti- 
mates reveals the presence of significant gas fractions (~ 
40% in average), which are compatible with a simple esti- 
mation from thei r star- format ion rate density assuming the 
iKennicuttl (| 19981 ) law. We observe typically lower gas frac- 
tions in the high mass objects. 

• We estimate the ionization parameter U for 8 objects 
where 032 and R23 line ratios are available, and derive high 
values with \og(U) ~ —2.5. The highest ionization values 
seem to correlate with the youngest stellar populations (< 
100 Myrs). 

• The gas-phase metallicites are calculated combining 
various line-ratios estimators, and we find a weaker evolu- 
tion in the mass-metallicity relation compared to estimates 
from bright LBGs observed in blank fields, with an offset 
reaching ~ 0.25 dex in the low stellar mass range (~ 10 9 
M©). This effect is seen both in the majority of individ- 
ual sources as well as in composite spectra created from the 
highest signal-to-noise or the low-mass objects. 

• Assuming that the evolution in the mass-metallicity re- 
lation is due to the increasing SFR at higher redshifts, we 
can reconcile our results with the existence of a fundamen- 
tal relation of mass, metallicity and SFR as proposed by 
iMannucci et al.l (|201Ch . The weaker evolution in the mass- 
metallicity relation in our sample is due to lower SFRs (com- 



pared to other luminous samples) for only slighly smaller 
masses. 

We foresee that the next HST programs on lensing 
clusters will continue to detect large number of magnified 
high redshift galaxies as multiple images, which would be 
ideal targets for deeper multi-object spectroscopic similar 
to the current work. This is a unique opportunity to extend 
the current sample to much lower stellar masses (typically 
10 8 Mq) and consequently provide more constraints on the 
mass-metallicity relation for a wide range of star-formation 
rates, metallicities and redshifts. 
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Figure 11. Variation between the measured meta llicity and 
the fu ndamental metallicity relation proposed bv lMannucci et all 
(2010), as a function of the stellar mass. Blue/black points are 
the 1.5 < z < 2.5 and 2.5 < z < 3.5 subsamples, respectively. The 
fundamental relation predicts the average metallicity of our lensed 
galaxies in both redshift ranges. 
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Table 1. The current sample of lensed galaxies. From left to right: astrometry, redshift and reference, magnification and reference. The 
sources below the line represent those drawn from the literature. 
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Run 


Date 


Seeing (") 


Photometric? 


A 


2005 October 13 


0.5 


Photometric 


B 


2006 July 24 


0.8 


Clear 


C 


2007 January 12 


1.0 


Clear 


D 


2007 May 3 


0.5-0.6 


Clear 


E 


2007 September 1 


0.5 


Photometric 


F 


2008 March 23 


0.4-0.5 


Photometric 


G 


2008 August 24 


0.5-0.9 


Clear 



Table 2. NIRSPEC observing runs and conditions 
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ID 
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rip 
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[IN 11 j 


A68-C1 
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7450 ± 90 
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MACS0744 
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1540 ± 40 


461 ± 23 
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3.042 
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200 ± 20 


850 ±29.6 






RXJ1347-11 
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1130 ±32 
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426 ± 123 


< 39 
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1.679 


E 


N5 








1140 ±34 


323 ± 19 


MACS0025 
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E 
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87 ±10 


< 66.7 


MACS0451 


2.013 


E 


N6 




1174 ±61 
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3170 ± 59 


729 ± 30 


MACS1423 


2.530 


F 


N6,N7 




123 ± 40 


120 ± 22 


313 ±39 


< 45.8 


RXJ1053 


2.576 


F 


N3,N6,N7 


510 ± 26 


624 ± 32 


1980 ±44.5 


1640 ± 42 


213 ±18 


A1689-Highz 


4.860 


F 


N7 


43 ±17 










A2218-Ebbels 


2.518 


F 


N3,N6,N7 


164 ± 13 


108 ±11 


170 ± 13 


273 ± 18 


< 21 


A2218-Flanking 


2.518 


F 


N6,N7 




180 ± 14 


402 ± 20 


444 ± 22 


< 14 


MACS0712 


2.646 


F 


N3,N6,N7 


< 32 


352 ± 20 


719 ± 40 


1070 ± 34 


225 ± 19 


C10949 


2.394 


F 


N3,N5,N7 


321 ± 20 


< 42.3 


819 ±30 


403 ± 20 


< 12 


A1835 


2.071 


F 


N5 




400 ± 20 


879 ± 46 


1700 ± 150 


111 ±20 


A773 


2.300 


F 


N6,N7 




268 ± 21 


108 ±11 


78 ±9 


126 ± 12 


A2218-Mult 


3.104 


G 


N4,N6 


1310 ±40 


871 ±31 


1920 ± 70 






A2218-Smm 


2.517 


G 


N6 




688 ± 29 


1030 ± 34 






A68-C4 


2.622 


G 


N6 




97 ±15 


< 22 






CL2244 


2.239 


ISAAC 


J,H,K 


370 ± 80 


< 42 


67 ±10 


463 ± 70 


< 20 



[Sn]( a ) 



( a ) Sum of the [Sn] A6717 and A6731 line fluxes. 

S/N Composite [Nil] A6584/Ha=0.16±0.03, [Sn] AA6717, 6731/Ha=0.11±0.04 

Low M* Composite [Nil] A6584/Ha;=0.13±0.10 

Table 3. Emission line measurements. Fluxes are given in units of 10 -19 ergs s _1 cm -2 . 
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ID 


z 


ri/2 
(kpc) 


A68-C1 


1.583 


1.09±0.18 


CEYE 


3.074 


1.75±0.21 


80CL0CK 


2.736 


1.47±0.38 


MACS0744 


2.209 


1.00±0.22 


Sextet 


3.042 


0.30±0.10 


RXJ1347-11 


1.773 


2.48±0.27 


A1689-Blob 


2.595 




C10024 


1.679 


10.0±1.2 


MACS0025 


2.378 




MACS0451 


2.013 


2.50±0.33 


MACS1423 


2.530 




RXJ1053 


2.576 


3.62±0.45 


A1689-Highz 


4.860 




A2218-Ebbels 


2.518 




A2218-Flanking 


2.518 


2.36±0.55 


MACS0712 


2.646 


0.75±0.23 


C10949 


2.394 


3.50±0.88 


A1835 


2.071 


1.52±0.07 


A773 


2.300 


0.38±0.05 


A2218-Mult 


3.104 


3.75±0.44 


A2218-Smm 


2.517 


1.86±0.40 


A68-C4 


2.622 




CL2244 


2.2399 




Composite 1 


2.37 ±0.34 




Composite 2 


2.52 ±0.10 





(km/s) 



(10 10 M ) 



M* 

(10 10 M ) 



SFR SFRcorr log(Z) E(B-V) 

(M yr- 1 ) (M yr' 1 ) 
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.34 
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,18 
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.05 
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8.42+° 


.48 
.13 


0. 


.13 




8.67+° 


.12 
.12 








8.64l° 


.19 

.38 







AC114-S2( b ) 

AC114-A2( b ) 

HORSESHOE^) 

CLONE^) 

J0900±2234( d ) 



1.867 

1.869 

2.38 

2.00 

2.032 



2.5 
2.9 



0.53±0.12 

2.36±0.67 

1.0 

2.2 

7.2 



32+ - 10 

U ' OZ -0.12 

52+ - 31 

u,oz -0.09 



0.6 



30. 
15. 
73. 
32. 
116±16 



Hjones et a l. (2010b) ( b ) |Lemoine-Busserolle et~aTl (120031 ) ( ^Hainline et al.l (l2009h ( d ^Bh 



8.78±0.20( e ) 
8.94±0.20 
8.49±0.16 
8.51±0.20 
8.12±0.19 
et al l hOld ) 



0.30 
0.40 
0.15 
0.24 
0.25 



( e )log(Z) value derived using the upper branch of the R23 diagram (see text for details). 

Table 4. Physical properties of the sample. From left to right: ID, redshift, magnification, half-light radius, measured velocity dispersion 
(corrected from instrumental resolution, see text for details), dynamical mass, stellar mass, SFR, metallicity, stellar extinction from the 
SED fitting. The SFR and stellar masses are corrected for the lensing magnification factor, and include the aperture corrections. Average 
redshift and stellar mass are given for the composite spectra described in the text. 




Figure Al. Central region around the BCG of the lensing cluster 
RXJ1053, showing the giant arc A visible as 4 images Al to A4. 
The red line outlines the critical curve at the redshift z = 2.576 
of the source A. 



APPENDIX A: RXJ1053 MASS MODEL 

In order to derive the magnification factor for the source at 
z = 2.576 in the lensing cluster RXJ1053, we constructed 
a parametric mass model of the central region of the clus- 
ter using the Lenstool software (I Julio et alJl2007h. Follow - 



ing similar strong- lensing works (e.g. Richard et al. I [20091 ), 
we assume the cluster mass distribu tion to follow a double 
Pseudo-Isothermal Elliptical fdPIE. llEKasdottir et afll2007h 
profile, and we add two central cluster members as lower- 
scale perturbations in the mass distribution. This model is 
constrained by the location of 3 images of the giant arc, 
clearly detected thanks to their morphology and symmetry 
on the V-band HST image fFig. lAl| . A fourth central image 
(A4) is predicted by the best fit model and identified on the 
same image. The best fit model has an rms a — 0.15" be- 
tween the predicted and observed positions of the 4 images. 
We report the best-fit parameters of the mass distribution 
in Table [ATI 
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Component x y e a r CO re r cu t 

(") (") (deg.) (kms- 1 ) (kpc) (kpc) 

Cluster [0] [0] 0.55±0.11 52.6±4.5 7051^ 60±23 1000.0 

BCG [0] [0] [0.146] [41.4] 705±150 [0.] 39±f 

Gall [2.1] [4.2] [0.16] [42.4] 202±87 [0.] 20±l° 

Table Al. Best fit parameters for the mass distribution reproducing the multiple components of the giant arc in RXJ1053. The dPIE 
parameters are given for each component: centre, ellipticity and orientation, velocity dispersion, core and cut radii. 



to 

Table A2. Multi-walength photometry (AB magnitudes) of the lensed sources in our sample. First column gives the UV (3 slope (see text for details). We use the notation BB'GVRII'JH l>0 
as a shortcut for the various HST broad-band filters. Filter names in italics refer to ground-based images. 

ID (3 F450W(B) F555W(G) F775W(F) F850LP F110W/J F160W/# K/Ks IRAC 3.6/im IRAC 4.5/im fe) 

F475W(B') F606W(V) F814W(I) 

F702W(R) © 
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